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We investigate Andreev reflection in two-dimensional hetero junctions formed by a superconductor 
in contact with a topological insulator ribbon either possessing or breaking time-reversal symmetry. 
Both classes of topological insulators exhibit perfect Andreev reflection, robust against disorder. 
This is assigned to topologically protected edge states. In the time-reversal symmetric case we show 
that doping one of the ribbon edges with magnetic impurities suppresses one Andreev channel, while 
no such suppression is seen in the broken symmetry situation. Based on this observation we suggest 
a transport experiment able to distinguish between the two types of topological insulators. 

PACS numbers: 



Topological insulators (TI) are a novel class of ma- 
terials, which have been intensively investigated in the 
past few years [H, 0). They come in two varieties, Z 2 
insulators, which preserve time-reversal symmetry, and 
Chern insulators, which do not. At present a variety of 
materials have been found to exhibit the time-reversal 
symmetric topological phase In contrast, although 

there are several theoretical proposals for constructing 
a Chern insulator none has been reported to date 

and the experimental quest continues. Also it is intrigu- 
ing the possibility of interfacing Z 2 insulators with su- 
perconductors (SCs). This interface has been predicted 
to host Majorana fermions, with possible applications in 
topological quantum computing [9[. 

A metal/superconductor interface may reflect an inci- 
dent electron from the metal as a positively charged hole 
with opposite spin, while a Cooper pair is formed in the 
superconductor. This electron- hole conversion is known 
as Andreev reflection [l(J and has long served as a use- 
ful probe for spin-polarized currents As topological 
insulators in general sustain edge states with a well de- 
fined spin, the Andreev reflection technique appears to 
be particularly suited to study edge states scattering in 
topological insulators. 

Recently there have been several interesting sugges- 
tions for turning the prototypical two-dimensional mate- 
rial graphene into a topological insulator. As far as the 
time-reversal symmetry broken case is concerned, Qiao et 
al. have shown, by means of first principles calculations, 
that Fe atoms adsorbed on graphene open up a band-gap 
in the bulk and yield a non-zero Chern number for the 
valence band Q. In another ingenious proposal Weeks 
et al. suggested to dope graphene with non-magnetic 
heavy adatoms in order to stabilize a time-reversal sym- 
metric topological insulating phase [l2| ■ In light of these 
promising developments, we have decided to study An- 
dreev reflection in two-dimensional topological insulator- 
superconductor junctions. 

We consider a two-dimensional topological insulator 
ribbon realized on a honeycomb lattice with zig-zag edge 



geometry, as shown in Fig. Q] The region to the right 
(SC region) is proximity coupled to a superconducting 
electrode, while the region on the left (TI region) is the 
topological insulator. The electron and hole spectra are 
described at the mean-field level by the Bogoliubov-de 
Gennes equation [l~3j . 
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where u and v are the wave-functions for electrons and 
holes, respectively. H is the single-particle Hamiltonian 
for the topological insulators, T is the time-reversal op- 
erator, A is the pairing potential and Ef is the Fermi 
level. In the left region (TI) the pairing potential is set 
to zero, i.e. there is no superconductivity. In the right 
region (SC) a finite constant pairing potential exists due 
to the proximity with a superconducting electrode. 

We use the Kane-Mele model [14j as single-particle 
Hamiltonian for the time reversal symmetric (Z 2 ) topo- 
logical insulator. This reads 
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Here the first term is just the nearest-neighbour hopping 
with strength t, where the spin indices of the creation, cj, 
and annihilation, a, operators have been omitted. The 
second term represents a staggered sub-lattice potential, 
i.e., the A type sub- lattice has an onsite energy A (£ = 
+1), while the B sub-lattice has onsite energy — A (£ = 
— 1). The last term describes second nearest-neighbour 
hopping with strength t 2 and it is purely imaginary (t 2 
is real and i = v— l). Furthermore z/y is equal to +1 for 
anti-clockwise hopping and to —1 for clockwise. Here a z 
is the z-component Pauli matrix describing the electron's 
spin. The last term can be thought as a mirror-symmetric 
spin-orbit interaction, since it couples the orbital motion 
of the electrons to their spins. 

For the time-reversal symmetry broken case we use a 
spinfull version of the Haldane model [15| . proposed by 
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Chen et al. [l6[ (from now on the spin-Haldane model, 
SH). The single particle hamiltonian reads 
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where the second term is the exchange field with strength 
7 , i.e. it represents Zeeman coupling. In addition to 
spin, also the orbital angular momentum of the electron, 
v>ij, is coupled to the exchange field. Following Chen et 
al. [16[ we approximate /3(j) w f3sgn(-f), and choose /3 to 



be negative. This parameter set describes a diamagnctic 
response to the magnetic field 7. Note that in this case 
the second nearest-neighbour hopping term has the same 
sign for both the spins, as opposed to that in Hkm ■ 




FIG. 1: Setup for calculating the two terminal transmission. 
Region SC is proximity coupled to a superconducting elec- 
trode while region TI is the topological insulator described by 
the two chosen single particle models. The rectangle marks 
the region at the TI/SC interface where disorder is intro- 
duced. 



We use the ballistic Landauer-Biittiker scheme [ItJ for 
calculating the transmission across the system. The self- 
energy matrix El (Sr) for the left-hand side (right-hand 
side) contact is obtained by using the electrodes' surface 
Green's function, g s . This is calculated iteratively from 
the following equation [l8j |. 
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where Hq is the Hamiltonian describing the electrodes 
unit cell and H\ is the coupling matrix between cells 
(note that in our tight-binding formulation the Hamil- 
tonian of the ribbon has a tri-diagonal form). The 
self-energy matrices are given by El,r = H^^g s H^ R , 
where H^ (Hn) is the coupling matrix between the scat- 
tering region and the left-hand side (right-hand side) 
contact. Then, the retarded Green's function, G r , for 
the scattering region of Hamiltonian -ffsc is obtained as 
G r = [{E + iO+)I - H sc - E L - Er]- 1 . The scattering 
region comprises the SC/TI interface and a portion of 
the electrodes. Finally the total transmission is simply 
T(E) = Tr(r L G r rRG rt ), where r LjR are the broadening 
matrices Tl.r = i[EL,R— E L R ]. Furthermore, the normal 



transmission coefficient from the n-th terminal to the ru- 
th one is obtained as T nma (E) = Tv{T n<J G r aa Y ma G r X) , 
while the Andrccv reflection coefficient is calculated as 



(5) 



where a = (7, \) and a = (J,, 7) are the spin indices. Thus 
R^u mS describes an incident electron from terminal n 
being reflected as an opposite spin hole into terminal m. 

We begin our analysis by calculating R A as a func- 
tion of energy, which is shown in figures [2]Ja) and (b) 
respectively for the Kane-Mele and spin-Haldane model. 
Henceforth, we set the nearest neighbor hopping t = 1 
and measure all the energies in units of t. Furthermore 
we fix the number of sites along the ribbon width to 
be n y = 18. The insets of Fig. [5] show the bandstruc- 
ture for the two models calculated in this strip geom- 
etry. In the bulk gap there exist gapless edge modes, 
a single pair on each edge. For the Hkm Hamiltonian, 
these are opposite spin Kramer's pairs, while for ffsH, 
there are two left-movers (one for each spin) at one edge 
and two right-movers at the other, exactly like integer 
Quantum Hall states. We find that for both these cases, 
the edge modes lead to a perfect Andreev reflection for 
electrons with energy smaller than the superconducting 
gap. In fact, in both cases the edge modes are perfectly 
Andreev reflected. Normal reflection, where an incident 
particle is reflected back without being converted into its 
antiparticle, is completely suppressed for the edge states 
as long as the Fermi energy lies in the bulk gap, as we 
have verified numerically These findings are consistent 
with recent theoretical and exnerimenta.l studies for time- 




FIG. 2: Andreev reflection coefficient for (a) Z2 and (b) Chern 
insulators showing perfect Andreev reflection for electron en- 
ergies smaller than the superconducting gap. The insets show 
the bandstructure for the two models solved in a ribbon ge- 
ometry. Here we choose £2 = 0.33, A/fe = 2.0, 7 = 0.20, 
f3 — —0.11 and A = 0.50. The fermi level Ef is set at zero. 



Next, we study whether such perfect Andreev reflec- 
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tion is robust to perturbations of the electronic structure 
at the SC/TI interface. To this goal we consider the 
effect of onsite disorder, which is introduced by adding 
a term of the form -ffdisordcr = Si u 'i c l c i to both Hkm 
and -ffsH- Here we choose the onsite energy, Wi, to be 
randomly distributed within the interval [— W/2,W/2]. 
Such disorder is introduced in a n x = 15 site-long re- 
gion near the SC/TI interface. From Fig. [3] it can be 
clearly seen that the Andreev reflection process is very 
robust against disorder. Even for moderately large dis- 
order (W ss 2.0t) R A remains perfectly quantized. This 
is attributed to the presence of the topological protected 
edge states, which are highly immune to impurities and 
disorder, and the situation is identical for both classes 
of topological insulators. For W > 2.0t fluctuations in 
R A begin to develop in the energy range where only edge 
states exist. As a result the Fano factor becomes non- 
zero. This signals a transition from ballistic to diffusive 
transport where backscattering is allowed and the edge 

□ 




FIG. 3: (color online) Effect of onsite disorder on the Andreev 
reflection coefficient for (a) Z 2 and (b) Chern insulators. The 
Andreev process is highly robust against onsite disorder and 
the crossover to dissipative transport occurs for W ~ 3.0t for 
Z2 insulator and W ~ 2.0t for Chern insulator. Here again 
we set t 2 = 0.33, \/t 2 = 2.0, 7 = 0.20, f3 = -0.11 and 
A = 0.50 and the Fermi level Ef is taken at zero. The curves 
are averaged over 640 random configurations. 



Finally we introduce magnetic impurities at one of the 
TI ribbon edges. The exchange coupling between the 
electron s pin and the impurities is incorporated into the 



model as [22 1 
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where S„ is n-th spin component of the magnetic impu- 
rity located at the edge site i, and J z and J11 are respec- 
tively the longitudinal and transverse exchange coupling. 
In the notation we have now explicitly re-introduced the 



spin index so that c\ a (ci a ) is the creation (annihilation) 
operator for an electron at site i with spin a. 

When one includes only the z-component of the ex- 
change coupling in the Kanc-Mele model, there is a shift 
of the up and down spin edge bands, by an amount pro- 
portional to the coupling J z but no band-gap opens in the 
edge state spectrum. For small values of J z , before the 
bulk band-gap closes, the system is in the time-reversal 
symmetry broken Quantum Spin Hall phase predicted by 
Yang et al. [H] • As a consequence, although we have lo- 
cally broken time-reversal symmetry, perfectly quantized 
Andreev reflection still occurs. This is because in the en- 
ergy range within the superconducting gap there is only 
the counter-propagating opposite spin channel available 
to normal reflection. 

In contrast, if we also include the transverse compo- 
nent of the exchange, i.e. we take Jy ^ 0, then a gap is 
opened at the edge where the magnetic impurities have 
been located. The destruction of the helical edge states 
at one of the two edges leads to a suppression of this 
channel, which results in the Andreev reflection coeffi- 
cient dropping from two to one, as shown in Fig. HJa). 
Such a reduction of R A from 2 to 1 is almost perfect 
except for some bulk contributions at energies approxi- 
mately equal to the superconducting gap. The situation 
for the spin-Haldanc model is different and the magnetic 
impurities produce no effect, regardless of the magnitude 
of J z and J|| . This is expected, since the topological pro- 
tection of the edge states for a Chern insulator continues 
to hold even in the absence of time-reversal symmetry. 
Consequently, no such suppression is observed and the 
Andreev reflection coefficient remains perfectly quantized 
to a value of two, as illustrated in Fig. HJb). 

Thus, we have shown that Andreev reflection mea- 
surements can characterize a topological insulator and 
distinguish it from a topologically trivial material. Per- 
fect Andreev reflection provides a signature for existence 
of topological edge states, although it is not unique to 
them. One can in fact envisage other systems, which dis- 
play similar perfect electron hole conversion, for instance 
a pair of ballistic nanowires. What is unique, though, is 
the tremendous immunity to disorder, which both types 
of topological insulators display. Furthermore, we also 
showed that the Z2 and Chern insulators respond differ- 
ently to presence of magnetic impurities. 

Based on this observation we propose a transport ex- 
periment to distinguish between the two types of topolog- 
ical insulators. The experiment involves placing magnetic 
impurities along one of the edges of the two-dimensional 
sample, for instance by using the tip of a Scanning Tun- 
neling Microscope. The impurities' spin will, in general, 
be aligned in arbitrary directions. The illumination with 
low-frequency polarized infrared light can however induce 
their alignment. This has been demonstrated, for in- 
stance, for Mn impurities in CdTe j24[. The infrared 
pulse imparts a momentum to align the impurity spins, 
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FIG. 4: Andreev reflection coefficient in presence of magnetic 
impurities located at the edge of a TI ribbon: (a) Zi insulator 
(b) Chern insulator. The suppression of one of the edge chan- 
nels in the time-reversal symmetric case R A produces a drop 
in R A from 2 to 1. Magnetic impurities have no effect on the 
Andreev reflection for a time-reversal symmetry broken insu- 
lator. Here we have chosen J, = = 0.50 and |S| = 2. The 
other parameters are same as before. 



which subsequently relax back to their random orienta- 
tions. The Andreev reflection coefficient R A can then 
be measured as a function of time, and this can be re- 
lated to the inclination angle 8 of the impurity spin S. 
The dependence of R A on 8 is shown in Fig. [5] As the 
spin rotates towards the z-direction R A returns back to 
the perfectly quantized value of two, the same as that in 
absence of impurities. For Chern insulator the Andreev 
reflection process is unaffected by magnetic impurities 
and thus to the exposure with polarized light. 
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FIG. 5: Andreev reflection coefficient in the presence of 
magnetic impurities for Z2 insulator as a function of the spin 
inclination angle 9 for various values of the exchange coupling. 



The relaxation of the spins leads to R 
unity. 



reverting back towards 



consists in applying an intense static magnetic field B 
perpendicular to the plane of the sample, and then to 
switch it off over a short time scale tg. The impurity 
spins will then relax back to their random configuration 
with a typical spin- relaxation time ts- During the time 
window comprised between ts and ts, measurements of 
R A should yield a behaviour similar to that shown in 
Fig. El We note that these possible approaches were al- 
ready outlined in Ref. relative to the observation of 
the Chern insulating phase in Mn doped HgTe quantum 
wells. 

In conclusion, we have investigated SC/TI heterojunc- 
tions and shown that they display perfect Andreev re- 
flection. The robustness of the topologically protected 
edge states lends this effect a large immunity against dis- 
order. We have then looked at magnetic impurities and 
shown that in the case of transverse exchange coupling 
the Andreev reflection coefficient of Z 2 topological insu- 
lators drops from two to one. This observation allowed 
us to propose a transport experiment able to distinguish 
between the two types of topological insulators. This 
consists in following the time evolution of the Andreev 
reflection coefficient of a device dusted with magnetic im- 
purities, which have been previously polarized. 
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A second possible route to spin-polarize the impurities 
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